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Abstract 

Artificial antiferromagnets and synthetic metamagnets have attracted much attention recently 
due to their potential for many different applications. Under some simplifying assumptions these 
systems can be modeled by thin Ising metamagnetic films. In this paper we study, using both the 
Wang/Landau scheme and importance sampling Monte Carlo simulations, the equilibrium prop- 
erties of these films. On the one hand we discuss the microcanonical density of states and its 
prominent features. On the other we analyze canonically various global and layer quantities. We 
obtain the phase diagram of thin Ising metamagnets as a function of temperature and external 
magnetic field. Whereas the phase diagram of the bulk system only exhibits one phase transition 
between the antiferromagnetic and paramagnetic phases, the phase diagram of thin Ising metamag- 
nets includes an additional intermediate phase where one of the surface layers has aligned itself with 
the direction of the applied magnetic field. This additional phase transition is discontinuous and 
ends in a critical end point. Consequently, it is possible to gradually go from the antiferromagnetic 
phase to the intermediate phase without passing through a phase transition. 

PACS numbers: 75.10.Hk,75.30.Kz,75.70.-i,75.40.Mg 
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I. INTRODUCTION 



Artificial antiferromagnets, formed of nanostructured superlattices that are coupled anti- 
ferromagnetically, have been the focus of many recent studies, due to their high potential for 
innovative technological applications. These possible applications range from high-density 
recording technology^ to spintronics devices^ and magnetic refrigeration.- Artificial antifer- 
romagnets are heterostructures composed of ferromagnetic layers that are coupled antifer- 
romagnetically via spacers. This structure yields a high level of control over both the intra- 
and interlayer interactions, allowing for a tailoring of the physical properties. Examples 
include [Co/Pt]/Ru, where ferromagnetic Co/Pt multilayers are periodically separated by 
Ru layers,^'^ or Co/Cr, where thin Co films are separated by the spacer Cr.- The magnetic 
moments of the ferromagnetic layers can thereby either be perpendicular (as it is the case 
for Co/Pt with perpendicular anisotropy) or parallel (as encountered in Co films) to the 
interface separating the ferromagnetic layers from the spacers. 

[Co/Pt]/Ru with strong perpendicular anisotropy, as well as the related systems like 
[Co/Pt] /NiO, Co/Ir, and Fe/Au, have been called synthetic metamagnets,- as they can 
exhibit a regime with an antiferromagnetic phase at low external magnetic fields and a 
paramagnetic phase at large fields. When changing the value of the magnetic field, plateaux 
show up in the magnetic hysteresis of these films, due to the fact that different layers reverse 
their magnetization at different values of the field.-^ In fact, for an even number of identical 
ferromagnetic layers in a thin film exactly two metamagnetic transitions are observed, as one 
of the outermost layers switches at a lower field than the internal layers.- The same sequence 
of phases should also be realized for in-plane magnetization and strong anisotropy.- 

In a phenomenological approach it is customary to replace in systems with strong per- 
pendicular anisotropy the ferromagnetic multilayers by a single ferromagnetic layer.-'-'' This 
naturally leads to a modeling of synthetic metamagnets by an Ising metamagnet, i.e. a lay- 
ered Ising model with ferromagnetic in-layer interactions and antiferromagetic interactions 
between adjacent layers. The advantage of using a layered Ising model is that this type 
of model is perfectly suited for a study of thermal properties through standard numerical 
methods. However, modeling ferromagnetic multilayers by single ferromagnetic layers also 
has its restrictions as it does not allow a theoretical description of the multidomain states, 
composed of metamagnetic stripe and bubble domains, that may form when applying a 
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magnetic field to an artificial antiferromagnet with strong anisotropy.-i^ 

Over the years Ising metamagnets have attracted quite some attention on their own.-"— 
Indeed, due to their simplicity, Ising metamagnets are ideal systems to study the properties 
of a tricritical point that separates a discontinuous transition between the antiferromagnetic 
and paramagnetic phases at high fields and low temperatures from a continuous transition 
between the same phases at low fields and high temperatures. In addition, the mean- 
field prediction of a decomposition of the tricritical point into a critical end point and a 
double critical end poinlr^ has led to systematic numerical investigations of this possible 
scenario. Whereas many studies have concluded that this decomposition does not take 
place in short-range models,— ii^i^Ii^^ the increase of the interlayer coordination, which brings 
the model closer to be of mean-field type, yields anomalies in the specific heat and the 
magnet izatioii^'i^ that resemble those measured in FeBr2.— 

The general behavior of Ising metamagnets is by now well understood. However, all men- 
tioned studies focused on bulk systems and no systematic studies of metamagnetic properties 
of thin Ising films have been done. The recent interest in synthetic metamagnets warrants 
an in-depth understanding of metamagnets in confined geometries, and we propose here a 
step in that direction by studying thin Ising metamagnetic films. 

An additional motivation for our work comes from the recent analysis^^ of non-equilibrium 
relaxation processes in Co/Cr superlattices. This study revealed that intriguing aging phe- 
nomena take place in layered antiferromagnets. In order to better understand these observa- 
tions a thorough study of the dynamical properties of related theoretical models is needed. 
However, before being able to study in depth the non-equilibrium properties of metamag- 
netic films, we found it necessary to first fully understand the equilibrium properties of 
these systems. Therefore we focus in this paper on the equilibrium phase diagram of thin 
metamagnetic films. The non-equilibrium properties of these systems will be discussed in a 
separate publication. 

Our paper is organized in the following way. After having introduced our model in Section 
II, we discuss in Section III its equilibrium properties. Using the Wang-Landau scheme, we 
determine the density of states (degeneracy) of our classical model and discuss its prominent 
features. This density of states is then used for a canonical analysis of the system where the 
focus is on the two phase transitions that are observed in thin film geometry when increasing 
the external magnetic field. As the Wang-Landau scheme is restricted to small systems, we 
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supplement our study by standard Monte Carlo simulations of larger systems. From these 
numerical data we derive the phase diagram of thin metamagnetic films and show that in 
thin films a new phase transition line shows up that separates two ordered phases. This 
transition, which is absent in bulk systems, is found to be discontinuous and to end in a 
critical end point. Finally, in Section IV we discuss our results and conclude. 

II. MODEL AND METHODS 

We consider a layered lattice model on a cubic lattice where every lattice point i is 
characterized by an Ising spin. Si = ±1. The interactions between nearest neighbor spins 
are ferromagnetic in the planes perpendicular to the z axis. These two-dimensional planes 
are coupled antiferromagnetically in the z direction. Adding an external magnetic field of 
strength B, the Hamiltonian of our system is then given by 



where we have introduced the internal energy E = —J^y ^i^j + JzYl ^i^j mag- 
netization M = ^2 Si. The sums over resp. are sums over nearest neighbor pairs 

i 

in the plane resp. along the z direction. The intralayer and interlayer coupling strengths 
are given by J^y > and > 0. 

In contrast to previous studies we focus on thin films which are realized by imposing 
free boundary conditions in the z direction, whereas in the x and y directions we have 
periodic boundary conditions. In accordance with the synthetic metamagnets we consider 
rather few layers, typically L = 8 or L = 10 (we restrict ourselves to even numbers of 
layers). For the case of periodic boundary conditions in all three directions, which has been 
studied extensively in the past, the system undergoes a metamagnetic transition between 
an antiferromagnetic phase and a paramagnetic phase when increasing the field strength. 
This transition is discontinuous for high fields and low temperatures and continuous for low 
fields and high temperatures. It is expected that an additional transition shows up when 
considering films, this transition being characterized by the alignment of the magnetization 
of one of the outermost planes with the external field.-"^ 

We study this system using different simulation techniques. In order to elucidate its 
static properties we compute the degeneracy (or density of states) Q{E, M) as a function of 




(1) 
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internal energy E and magnetization M using the Wang-Landau scheme.— The degeneracy 
can then be used for a standard canonical analysis, with the partition sum as a function of 
temperature T and magnetic field B being given by (we choose units for which ks = 1) 

Z{T, B) = J2Y1 M)e-^/'^+^^/'^ . (2) 

E M 

All global quantities of interest (like the mean magnetization, the mean energy, the suscep- 
tibility, and the specific heat) then follow from derivatives of Z with respect to T and B. 
It is well known that the Wang-Landau scheme usually yields a very good estimate for the 
degeneracy, allowing for a detailed canonical analysis. However, only rather small systems 
can be studied in this way. We therefore supplement our study with traditional importance 
sampling simulations of larger system sizes, using the Metropolis algorithm. We thereby 
study thin films with layers that contain L x L spins, with L = 8, 16, 32, and 64. Relevant 
quantities are computed for a large number of temperatures and field strengths; in order to 
completely cover the phase diagram, the increment between successive T and B values is 
typically 0.01. 

III. EQUILIBRIUM PROPERTIES 

We focus in the following on the thermal equilibrium properties of thin Ising ferromag- 
nets. The statistical physics treatment of our classical spin models allows for an in-depth 
discussion of its properties as a function of all relevant parameters, which in our case are 
the temperature and the strength of the magnetic field. 

A. The density of states 

As already mentioned we compute for our smaller systems the density of states Q{E, M), 
i.e. the number of microscopic configurations that have the same internal energy E and 
the same magnetization M. Having this two-dimensional histogram at our disposal, we can 
then compute all relevant global thermodynamic quantities from the partition sum ([2]) and 
its derivatives by simply inserting the numerical values for T and B. 

Before doing that we briefly discuss the density of states itself. Fig. [1] shows the micro- 
canonical entropy S{E, M) = lnQ{E, M) as a function of the energy density e = E/N and 



5 



the magnetization density m = M/N (here is the total number of spins in our systems) 
for three different cases. Whereas in (a) we consider a system with periodic boundary con- 
ditions composed ofA^ = 8x8x8 spins, with J^y = J^, in (b) and (c) we show two systems 
ofA^ = 8x8x8 sites with free boundary condition in z direction, the different systems 
having different relative strengths of the interactions: J^y = Jz in (b) and J^y = 2Jz in (c). 




-3 -2 -1 1 2 3-3 -2 -1 1 2 3-3 -2 -1 1 2 3 

E/N 



FIG. 1: Microcanonical entropy S{e, m) as a function of energy density e = E/N and magnetization 
density m = M/N for systems containing = 8 x 8 x 8 spins, (a) Periodic boundary conditions 
in all three directions and Jxy = Jz, (b) free surfaces in z direction with Jxy = Jz, and (c) free 
surfaces in z direction with J^y = 2Jz- The increment between two successive contour lines is 20. 
The supports of the entropy surfaces are indicated by the larger points. 

In a microcanonical analysis one infers the physical property of a system from a direct 
study of the microcanonical entropy ^(e, m).—-^ Most investigations of this type focused on 
spin models with ferromagnetic nearest neighbor interactions, as for example the standard 
Ising or Potts models,—"— or on polymer models. Due to its complicated interactions, 
the microcanonical entropy of the Ising metamagnet is more complex as, for example, that 
of the standard nearest neighbor Ising model,— see Fig. [1] Still some of its features and 
properties are readily understood. 

Firstly, there are obvious properties that are independent of the boundary conditions and 
of the interaction ratio Jz/ Jxy The ground state is two-fold degenerate and has magneti- 
zation zero (recall that we only study systems with even numbers of layers), with the fully 
ordered layers pointing alternatively in up or down direction. Obviously, a given ground 
state can be changed into the other ground state by multiplying all the spins by —1. This 
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symmetry of the internal energy also shows up in the symmetry of the entropy surfaces with 
respect to the M = line. Interestingly, the entropy surfaces all display prominent kinks at 
magnetizations M/N = ±2/i/L, with /i = 0, ■ ■ ■ , L/2, corresponding to configurations with 
fully ordered planes where fi planes have been flipped with respect to the ground state. 

Changing the interaction ratio Jz/Jxy, see Fig. [T]d and dt, mainly changes the range of 
accessible energies. As a result, smaller values of Jz/Jxy yield a compressed, but otherwise 
unchanged, entropy surface. 

At a first look, changing the boundary condition in z direction also seems to only have 
minor effects on s{e,m). However, one notes that the energy difference between the ground 
state and the state with one flipped layer, with M/N = ±2/L = ±1/4 in Fig. [H is smaller for 
the free boundaries in z direction than for the periodic boundaries, compare Fig. [T^ and[DD. 
As we already argued in the introduction (and as we will see later in the canonical analysis), 
our metamagnetic films should display two metamagnetic transitions, a first transition at low 
magnetic fields, where only one surface layer is flipped, and a second transition at higher 
fields where the remaining layers pointing opposite to the magnetic field are flipped. In 
contrast, a system with periodic boundary conditions in all three directions is known to 
undergo a single metamagnetic transition.— With this in mind, it seems surprising that very 
similarly looking entropy surfaces should yield these different phase transition sequences. 

In fact, it is the decrease of the energy difference between the ground state and the states 
with fully ordered planes and magnetization M/N = ±2/L that ultimately is responsible 
for the emergence of this additional transition. We show that in Fig. [2] where we plot for 
systems with J^ = Jxy the quantity {E — BM — TS)/N as a function of M/N, with T = 1 
and different values of B. In the canonical ensemble the corresponding quantity is of course 
the Helmholtz free energy density which is minimal for the stable state. But even without 
making a canonical average one can read off the stable phase from the "microcanonical" 
quantity plotted in Fig. O For the system with periodic boundary conditions, see Fig. 
Ell, our quantity is minimal for M = and low fields. When the field strength exceeds 
B = 2, a metamagnetic phase transition takes place such that the stable phase is now the 
paramagnetic phase with M/N = 1. Conversely, for the thin film geometry, see Fig. [2}3, 
a first transition to a phase with one flipped layer and magnetization M/N = 2/L = 1/4 
shows up at i? = 1, as can be seen when studying the global minimum of E — BM — TS, 
followed by a second discontinuous transition at B = 2. In this way, the sequence of phase 
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transitions can indeed be unraveled in a microcanonical analysis. 




FIG. 2: (Color online) "Microcanonical" quantity {E - BM - TS)/N as a function of M/N. 
The canonical correspondence of this quantity is the Helmholtz free energy density. The global 
minimum of the plotted quantity allows to understand the sequence of phases that show up when 
changing the strength of the magnetic field. The systems contain = 8 x 8 x 8 spins, with (a) 
periodic boundary conditions and (b) free boundary conditions in z direction. The insets show the 
relevant parts of the microcanonical entropy, see main text. 



B. Thermal quantities 

Even so the microcanonical analysis allows to determine the sequence of phases, an in- 
depth study of the thermal properties of the metamagnetic films warrants a canonical anal- 
ysis of the standard quantities as a function of temperature and magnetic field. In the 
following, we discuss small systems for which the density of states can be obtained numer- 
ically, as discussed in the previous subsection. This density of states M) = e^^^''^^^ 
is then used for the computation of the partition sum and related quantities. Thus for a 



8 



quantity Q{E, M) the thermal average at temperature T and field strength B is 

E M 

where the partition function Z{T, B) is given by Eq. ([2]). In order to elucidate the thermal 
properties of thin metamagnets we studied the average magnetization (M), the average 
energy (£'), the magnetic susceptibility x = [(^^) ~ (^)^]) the specific heat C = 

We supplement this canonical analysis of small systems based on the density of states 
by standard canonical simulations of larger systems, with the same number of layers but 
different number of spins in the layers. This allows us to assess the finite-size effects and 
to extrapolate to films with a large number of spins in every layer. Standard importance 
sampling Monte Carlo simulations often yield results of lesser quality than the canonical 
analysis based on the Wang/Landau scheme. On the other hand larger systems can easily 
be simulated. In addition, we also can look at additional quantities which are not accessible 
when starting from the degeneracy as a function of total magnetization and total energy. 
Thus we will also study the layer susceptibility 

X{z) = ^ [{M\z)) - {M{z)f] (4) 

where M{z) is the magnetization of layer z. Of special interest is of course the susceptibility 
of the layer z = 1 that fiips in the direction of the magnetic field when crossing the low 
field phase transition line. Similarly, we also study the layer specific heat 

c{z) = [{n\z)) - {n{z)r] (5) 

layer z 

where T-Liz) = ~Jxy Yl ^i^j ^ ^ Yl '^'j is the in-layer contribution to the energy in layer 

Z. 

Fig. [3] highlights the expected difference in the field dependence of the magnetization 
between thin films and bulk systems (i.e. systems with periodic boundary conditions). For 
thin films a first transition takes place for B ^ J^, where one of the outer layers completely 
fiips, followed by a second transition at i? ^ 27^, where the remaining layers align with the 
magnetic field. At T = these transitions take place exactly at B = and B = 2J^. For 
larger T than that used in Fig. [3l the total magnetization after the fiipping of the outer 
layer is slightly lower than N/4, due to thermal fiuctuations. 
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FIG. 3: (Color online) Average magnetization as a function of magnetic field strength B at tem- 
perature T = 1. The magnetization is shown for thin metamagnetic films with different interaction 
ratios as well as for a sample with periodic boundary conditions (pbc) and = 1. The three- 
dimensional samples contain = 8^ spins, with the interlayer interaction having the strength 
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We discuss the T and B dependence of our thermal quantities in Fig. |4] for a thin 
film of eight layers with = Jxy = 1- The two phase transitions are readily seen in the 
changes of the magnetization density (plateaux in Fig. Ht) and the susceptibility (lines of 
maxima in Fig. HJi). One notes that for small T the change in the magnetization is very 
abrupt, pointing to a discontinuous transition, whereas for larger T this change is much 
more gradual, in agreement with a continuous transition. This change of the nature of the 
transition also shows up in the susceptibility. Two lines of maxima, one being hardly visible 
on the scale of the figure, are also observed in the specific heat, see Fig. |3)d. 

Based on the positions of the maxima in the response functions, see Figures O U\ and[8l we 
obtain the phase diagram for a thin metamagnetic film composed of 8 layers shown in Fig. 
We first note the expected presence of three different phases: the paramagnetic phase at high 
temperatures and high fields, the antiferromagnetic phase at low temperatures and low fields. 
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FIG. 4: (Color online) (a) Energy density, (b) specific heat, (c) magnetization density, and (d) 
susceptibility as a function of T and B for a thin film composed of = 8^ spins, with Jz = Jxy = 1- 

and a phase which has a non-zero magnetization at intermediate fields and low temperatures. 
Similar to what is observed in the phase diagram of the bulk system, the transition between 
the ordered and paramagnetic phases is discontinuous at low temperatures and continuous 
at high temperatures, with a tricritical point separating these two regimes. We locate this 
point at T = 2.20(1) and B = 1.96(1). At low temperatures the transition between the 
antiferromagnetic phase and the intermediate phase is also discontinuous. This transition, 
however, does not extend all the way up to the phase transition line separating the ordered 
phase from the paramagnetic phase, but instead ends at a critical end point located at 
T = 2.50(5) and B = 0.98(1), see below. 

Increasing the thickness of the film leaves the phase diagram qualitatively unchanged. A 
slight shift in the phase transition lines, especially at high temperatures and low fields, is ob- 
served when changing the number of layers in the film. This shift of the critical temperature 
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FIG. 5: (Color online) Phase diagram for the thin Ising metamagnet film with eight layers. The 
different transition lines result from the maxima of the (global and layer) susceptibility and specific 
heat. Black lines indicate a discontinuous transition, whereas the continuous transition is shown 
as a grey (cyan online) line. We observe a tricritical point (filled point) as well as a critical end 
point (open point). 

of a film as a function of thickness is of course expected and has been studied extensively, 
both theoretically and experimentally, in the absence of external magnetic fields (see, for 



In Figures |6l U\ and |8] we have a closer look at the different response functions used 
for the construction of the phase diagram. At low temperatures, see Fig. [6^, the peaks 
in the susceptibility as a function of the magnetic field strength are very pronounced and 
sharp, indicating the discontinuous character of the transitions. The form of the peaks 
change around T = 2.2, as here the order of the transition changes from discontinuous to 
continuous. Above T = 2.2, see Fig. the height of the peaks shows the expected size 
dependence of a continuous transition. Fig. [6li shows the total susceptibility at the rather 
high temperature of T = 3.5. Increasing the system size reveals the emergence of a critical 
peak at B ^ 1.3. This peak, which coincides which the transition to the paramagnetic 
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FIG. 6: (Color online) Susceptibility of thin films with 8 layers as a function of magnetic field 
strength for various temperatures. (a,b) Data derived from the density of states of a small system 
with 8x8 spins per layer. At low temperatures the phase transition is discontinuous, as revealed 
by the characteristic form of the peaks. Around T ~ 2.2 the form of the peaks changes, indicating 
that the transition has become continuous, (c) Susceptibility for T = 2.5 obtained for systems with 
different numbers of spins in the layers. Whereas the height of the peaks changes with system size, 
the location of the peaks does not show a dependence on the sizes of the layers. Symbols indicate 
data obtained using the Metropolis algorithm, (d) The same as in (c), but now for T = 3.5. For 
larger systems a critical peak, which indicates the transition to the paramagnetic phase, emerges 
that is not visible for the smallest systems. The peak at i? ~ 1.6 is a non-critical peak. Error bars 
are comparable to the sizes of the symbols. 

phase, only shows up as a shoulder in the smaller systems. The additional peak at i? ~ 1.6 
is a non-critical one and is similar to that observed in the bulk system. 

Another way to monitor the phase transition lines is through a study of the specific heat 
as a function of temperature. As an example Fig. [7^ shows the specific heat for vanishing 
magnetic field. Only one critical peak is observed that results from the phase transition 
between the ordered low temperature phase and the disordered high temperature phase. In 
Fig. [7)d we show the field dependence of the maximum of the specific heat along the phase 
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FIG. 7: (Color online) (a) Specific heat of an eight-layer system as a function of T for = and 
various sizes of the layers. The peaks indicate the transition between the ordered and paramagnetic 
phases, (b) The maximum height of the specific heat for a fixed value of the magnetic field strength 
B. This maximum value is achieved at the transition to the paramagnetic phase. The strong peak 
around B ~ 1.96 reveals the change of the order of the transition when crossing the tricritical 
point. 

transition line to the paramagnetic phase. For low fields this transition is continuous and 
the specific heat height shows the expected finite size scaling of an ordinary critical point. 
For a fixed system size the specific heat exhibits a strong peak around B ^ 1.96 which is 
due to the change of the order of the phase transition when crossing the tricritical point. 

The standard way to locate the tricritical point is to monitor the hysteris observed when 
crossing the phase transition line at low fields^i2^, as the hysteresis loop vanishes when 
approaching the tricritical point. We found it useful to monitor in addition the location 
of a non-critical high field peak observed for temperatures larger than the temperature of 
the tricritical point. This non-critical peak is also observed in metamagnetic bulk systems, 
both in simulationa^S and in experiments^. The merging of this non-critical line with the 
phase transition line separating the paramagnetic phase from the ordered phases allows to 
reliably estimate the location of the tricritical point. Another estimate can be obtained by 
monitoring the strong increase of the peak heights of response functions when approaching 
the tricritical point, as shown in Fig. [7)d. Taking all this into account, we estimate the 
tricritical point of our thin film to be at T = 2.20(1) and B = 1.96(1). 



14 



T— I— I— |— I— I— I— I— |— I— I— I— r 




0.5 1.0 1.5 2.5 3.0 3.5 

B T 

FIG. 8: (Color online) (a) Surface susceptibility ([!]) at T = 3 and (b) surface specific heat 
at B = 0.5 for the surface layer that aligns with the magnetic field. Data obtained for different 
system sizes are shown. These maxima are non-critical and size independent, indicating that the 
reversal of the surface magnetization is gradual for these values of the system parameters. 

An interesting property of the phase diagram of thin metamagnetic films emerges when 
studying response functions at high temperatures and low fields. Whereas for low temper- 
atures the response functions display a characteristic peak with a height that depends on 
the system size, for higher temperatures the response functions do no longer display any 
dependence on the size of the system. As shown in Fig. [8] for the surface response functions, 
system size dependent peaks are no longer encountered in the ordered phase for temper- 
atures T > 2.50 and fields B > 0.98. This excludes the existence of a line of continuous 
phase transitions at high temperatures and low fields, indicating that the line of discontinu- 
ous transitions ends in a critical end point. Therefore, instead of crossing a phase transition 
line, with its characteristic singularities, one can go from the antiferromagnetic phase at low 
temperatures and fields to the intermediate phase with one flipped layer in a smooth and 
gradual way. This is in complete analogy to the well known behavior of water where one 
can go in a similar way from vapor to liquid without a phase transition. Based on our data 
we locate this critical end point at T = 2.50(5) and B = 0.98(1). 
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IV. DISCUSSION AND CONCLUSION 



The numerous possible applications have recently yielded a strong interest in artificial 
antiferromagnets and synthetic metamagnets. Remarkably, the layered structure of these 
materials allows for a high level of control of their physical properties through the fine-tuning 
of the strengths of both the inter- and intralayer interactions. 

In order to get a better understanding of some of these properties we presented a study of 
the equilibrium properties of thin Ising metamagnetic films. Indeed, neglecting the internal 
structure of ferromagnetic multilayers, one can model artificial antiferromagnets with strong 
perpendicular anisotropy by a layered Ising antiferromagnet. If one notes in addition that 
the synthetic metamagnets are usually composed of only a few repetition of the superlattice 
structure, one is then naturally led to consider Ising metamagnets in thin film geometry. 

Our work allowed us to study in detail the phase diagram of these systems as a function 
of temperature and of the strength of the applied external magnetic field. Earlier studies 
of related phenomenological models revealed^i^ that for an even number of layers a layered 
antiferromagnet has three different T = states, depending on the strength of the magnetic 
field. The antiferromagnetic structure, stable at low fields, goes over into a different ordered 
structure at intermediate field strengths. In this intermediate phase one of the surface layers 
aligns with the direction of the magnetic field, the other layers remaining unchanged. At 
larger fields the intermediate phase is replaced by the paramagnetic phase where all layers 
are parallel to the magnetic field. This intermediate phase is absent in the corresponding 
bulk system and is therefore characteristic for thin metamagnetic films. 

Using the Wang-Landau algorithm and importance sampling Monte Carlo simulations 
we thoroughly studied the equilibrium properties of our system. The Wang-Landau method 
yields the microcanonical entropy (or, alternatively, the degeneracy) as a function of magne- 
tization and internal energy. The microcanonical entropy surface is rather complicated and 
reveals interesting features that we discussed in some detail. Thus we were able to relate 
the appearance of the intermediate phase to a subtle change in the microcanonical entropy 
when changing the boundary condition. 

Our canonical analysis of different thermal quantities, as for example the global and layer 
magnetization, the global and layer susceptibility as well as the global and layer specific heat, 
allows us to determine the temperature — magnetic field phase diagram shown in Fig. |5l 
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which is the main result of our study. Interestingly, the discontinuous phase transition 
between the antiferromagnetic and the intermediate phase ends in a critical end point. 
It is therefore possible to pass from one phase to the other without undergoing a phase 
transition. Experimental studies on systems with strong perpendicular anisotropy should 
be able to verify this intriguing feature of thin Ising metamagnets. 

As already mentioned in the Introduction one of the motivations for our study was the 
recent investigation^^ of non-equilibrium relaxation and aging phenomena in Co/Cr su- 
perlattices. This study revealed intriguing non-equilibrium features that warrant a better 
understanding. As a first step in that direction we started a study of the non-equilibrium 
properties of thin Ising metamagnets, but soon realized that a more complete understand- 
ing of the equilibrium properties of these systems is needed before being able to develop 
better insights into the more complex situation of relaxation far from equilibrium. With 
the new knowledge of the equilibrium properties reported in this paper we are now well 
prepared to better understand these complicated non-equilibrium processes. The results of 
this investigation will be the subject of a forthcoming paper. 
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